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A  novel  and  unique  design  of  thermoelectric  generators,  in  which  a  heat  source  is  combined  with 
thermoelectric  elements,  is  proposed.  By  placing  heat-generating  radioactive  isotopes  inside  the  thermo¬ 
electric  elements,  the  heat  transfer  limitation  between  the  generator  and  the  heat  source  can  be  elimi¬ 
nated,  ensuring  simplicity.  The  inner  electrode  is  sandwiched  between  identical  thermoelectric 
elements,  which  naturally  allows  the  inner  core  to  act  as  the  hot  side.  Analysis  shows  that  conversion  effi¬ 
ciency  and  power  density  increase  as  the  heat  density  inside  the  thermoelectric  elements  increases  and 
as  the  thermoelectric  performance  of  the  material  improves.  The  theoretical  maximum  efficiency  is 
shown  to  be  50%.  However,  realistic  performance  under  practical  constraint  is  much  worse.  In  realistic 
cases,  the  efficiency  would  be  about  3%  at  best.  The  power  density  of  the  proposed  design  exhibits  a  much 
more  reasonable  value  as  high  as  3000  W/m2.  Although  the  efficiency  is  low,  the  simplicity  of  the 
proposed  design  combined  with  its  reasonable  power  density  may  result  in  some,  albeit  limited,  potential 
applications.  Further  investigation  must  be  performed  in  order  to  realize  such  potential. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 


Thermoelectric  generators,  which  utilize  thermoelectric  materi¬ 
als  to  convert  heat  to  electricity,  provide  an  alternative  option  to 
conventional  heat  engines  in  the  field  of  power  generation,  most 
notably  in  compact  waste  heat  recovery  systems  or  in  generic 
power  supplies  at  remote,  isolated  locations.  Typically,  the  maxi¬ 
mum  efficiency  of  thermoelectric  energy  conversion  is  presented 
in  terms  of  the  temperature  of  each  heat  reservoir  and  the  thermo¬ 
electric  figure  of  merit  ZT  =  S2T/(pk),  where  S  is  the  Seebeck  coef¬ 
ficient,  p  is  the  electric  resistivity,  k  is  the  thermal  conductivity, 
and  T  is  the  temperature,  leading  to  the  following  traditional 
expression  for  the  maximum  thermoelectric  conversion  efficiency 

too)  m; 


a) 


where  ZT  is  typically  evaluated  at  the  mean  temperature  (Tm)  of  the 
hot-side  temperature  Th  and  the  cold-side  temperature  Tc.  From  the 
equation,  it  is  clear  that  the  increase  of  ZT  leads  to  the  increase  of 
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the  energy  conversion  efficiency,  and  significant  efforts  have  been 
exerted  on  creating  materials  with  high  ZT  values.  Nanowires  [2], 
thin  films  [3],  and  materials  with  grain  boundaries  [4]  have  also 
been  investigated,  showing  various  degrees  of  success  in  improving 
ZT. 

However,  conventional  thermoelectric  generators  do  not  per¬ 
form  as  well  as  promised.  Many  factors  are  contributing  to  the  per¬ 
formance  gap  between  the  theoretical  estimate  of  thermoelectric 
materials  and  the  actual  performance  of  conventional  thermoelec¬ 
tric  generators.  One  of  the  most  significant  contributors  is  the  lim¬ 
ited  heat  transfer  between  the  hot-  and/or  cold-side  thermal 
reservoirs  and  the  generator  itself  [5],  Large  temperature  differ¬ 
ence  typically  occurs  between  the  heat  source  (or  sink)  and  the 
generator  due  to  a  poor  heat  transfer  coefficient.  The  actual 
temperature  difference  across  thermoelectric  elements  made  of 
thermoelectric  materials  becomes  small  as  a  result.  With  a  small 
actual  temperature  difference,  the  performance  of  thermoelectric 
materials,  which  is  relatively  poor  already,  suffers  even  more  in 
real-world  applications. 

In  this  article,  we  study  the  feasibility  of  a  self-heating  thermo¬ 
electric  generator,  in  which  thermoelectric  materials  create  their 
own  heat.  Since  the  thermoelectric  elements  heat  themselves,  there 
is  no  need  for  engineering  a  good  thermal  interface  between  the 
hot-side  thermal  reservoir  and  the  generator,  which  eliminates 
the  difficulty  of  achieving  good  heat  transfer,  at  least  on  one  side. 


http://dx.doi.org/10.1016/j.enconman.2014.02.061 
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Self-heating  thermoelectric  materials  are  not  being  utilized  yet,  but 
there  does  exist  certain  prospect  for  such  development.  One  poten¬ 
tial  way  to  develop  such  materials  is  to  incorporate  radioactive  iso¬ 
topes  in  traditional  thermoelectric  materials.  For  example,  a  known 
thermoelectric  oxide  (SrTi03)  [6,7]  can  be  altered  into  a  self-heating 
thermoelectric  oxide  by  replacing  regular  atoms  (Sr)  with  radioac¬ 
tive  isotopes  (90Sr).  Another  method  is  to  fill  radioactive  isotopes 
into  a  skutterudite,  which  is  one  of  the  most  promising  thermoelec¬ 
tric  materials  [8],  Radioactive  isotopes  were  previously  used  as  a 
heat  source  in  the  context  of  thermoelectric  generation  in  space 
missions  [9,10]  and  in  isolated  monitoring  sites  [11].  In  particular, 
the  use  of  90Sr  in  the  form  of  SrTi03  as  a  heat  source  was  already  re¬ 
ported  [12],  but  combining  radioactive  isotopes  with  thermoelectric 
elements  is  a  possibility  that  has  been  hardly  explored  before. 

The  paper  is  organized  as  follows.  In  Section  2,  first,  we  describe 
the  proposed  design  of  a  thermoelectric  generator  using  self-heat¬ 
ing  thermoelectric  materials.  Then,  we  continue  to  discuss  our 
method  of  analysis  for  evaluating  the  performance  of  the  genera¬ 
tor.  In  Section  3,  the  method  described  in  Section  2  is  applied  while 
changing  the  design  parameters  within  realistic  ranges  in  order  to 
evaluate  the  performance  indices  of  the  proposed  design.  In  Sec¬ 
tion  4,  conclusions  are  drawn  by  summarizing  and  discussing  the 
results. 

2.  Analysis  of  the  proposed  design 

Our  proposed  design  is  presented  in  Fig.  1.  Unlike  conventional 
thermoelectric  modules,  there  is  no  external  heat  source  other 
than  the  self-heating  thermoelectric  elements  inside  the  generator. 
Hence,  there  is  no  distinction  between  the  hot  side  and  the  cold 
side  specified  by  external  thermal  reservoirs.  Only  with  appropri¬ 
ate  thermal  boundary  conditions  imposed,  can  a  proper  tempera¬ 
ture  difference  across  a  thermoelectric  element  be  developed. 
The  use  of  a  sandwich-like  shape,  where  the  inner  electrode  is  lo¬ 
cated  between  two  identical  thermoelectric  elements,  is  used  to 
achieve  such  boundary  conditions  in  a  straightforward  fashion. 
The  inner  boundary  becomes  naturally  adiabatic,  which  renders 
the  inner  core  as  the  hot  side. 

For  simplicity,  we  only  consider  one  single  unipolar  element  in 
our  analysis.  Since  the  entire  generator  is  constructed  by  alternat¬ 
ing  concatenation  of  two  different  types  of  unipolar  elements,  the 
performance  of  the  entire  generator  can  be  analyzed  by  summing 
the  power  output  of  all  elements.  The  governing  equation  describ¬ 
ing  the  temperature  distribution  (T)  in  one  single  unipolar  element 
is  given  as  follows: 

4(k|)+/>/  +  «.0.  (2) 


x  is  the  variable  for  the  outward  longitudinal  coordinate,  k  and 
p  are  the  thermal  conductivity  and  the  electrical  resistivity  of  the 
thermoelectric  element,  respectively,  j  is  the  current  density,  and 
q  is  the  amount  of  heat  generated  per  unit  volume,  which  is  uni¬ 
form  and  constant  between  x  =  0  and  L.  Because  the  purpose  of 
the  analysis  is  a  preliminary  evaluation  of  feasibility,  only  minimal 
functional  features  are  included  in  Eq.  (2).  For  instance,  the  Thom¬ 
son  effect  is  excluded,  although  it  may  affect  actual  performance 
[13].  Another  effect  neglected  here  is  that  of  interfacial  Ohmic 
heating,  which  is  an  important  mechanism  causing  losses  in  tradi¬ 
tional  thermoelectric  modules,  where  the  leg  length  of  each  ther¬ 
moelectric  element  is  relatively  small,  like  several  millimeters.  In 
the  present  design,  however,  the  leg  length  must  be  much  larger, 
since  each  leg  must  act  as  its  own  heat  source  as  well.  As  shown 
in  Section  3,  the  leg  length  of  each  thermoelectric  element  may 
reach  more  than  5  cm.  Using  the  value  of  resistivity  in  Table  1, 
the  value  of  the  electrical  resistance  of  one  leg  per  unit  area  is  esti¬ 
mated  to  be  15  x  10  3  Q  cm2.  It  was  reported  that  the  interfacial 
resistance  of  a  typical  thermoelectric  module  is  on  the  order  of 
10  6  Q  cm2  [14],  Thus,  it  is  clear  that  the  interfacial  Ohmic  heating 
may  be  neglected  without  too  much  loss  of  generality. 

The  corresponding  boundary  conditions  are  given  by: 

jST-k^  =  0  at  x=0,  (3) 

and 

T  =  TC  at  x  =  L.  (4) 

Here,  S  is  the  Seebeck  coefficient  of  the  thermoelectric  element, 
L  is  the  coordinate  of  the  outer  boundary,  and  Tc  is  the  temperature 
of  the  cold-side  end.  The  first  term  on  the  LHS  of  Eq.  (3)  represents 
the  Peltier  heat.  From  Eqs.  (2)-(4),  dimensionless  equations  are 
constructed  by  defining  0  =  T/V  and  f  =  x/L.  T*  is  a  characteristic 
temperature,  which  is  given  here  as  T*  =  qL2/k.  The  set  of  corre¬ 
sponding  dimensionless  equations  are: 


^+v  +  l=0, 
dc2 

(5) 

Ve~%  =  0  at  i  =  °> 

(6) 

and 

e=ec  at  |=1, 

(7) 

where  0C  =  Tc/T*,  /( =jSL/k,  and  v  =  pj2/q.  The  solution  to  Eqs.  (5)- 
(7)  is  given  as  follows: 

(1  +  /2<f)(l  +  v  +  29c)  -  (1  +  /r)(l  +  v){2 
m  - '  (8) 


Fig.  1.  Schematic 


;  thermoelectric  elements. 
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Material  properties 


S(pV/K) 

-200 


k(W  m-'  IT1) 
3.5 


q  (W/cm3) 
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Once  the  solution  is  obtained,  the  temperature  distribution  can  be 
recovered  by  converting  the  result  into  a  dimensional  form.  Note 
that,  once  the  material  properties  (fc,  q.  S,  and  p)  and  the  geometry 
of  the  generator  (L)  are  specified,  p  and  v,  which  have  dependency 
on  the  operating  condition  (j),  are  not  independent  of  each  other, 
since 


(9) 


In  order  to  compute  the  power  output  per  axial  unit  length,  one 
must  specify  the  condition  of  the  external  load.  The  external  load 
is  modeled  as  a  simple  electrical  resistance  per  unit  area  of  the  gen¬ 
erator  (Rl),  which  can  be  adjusted  to  ensure  the  best  performance. 
We  define  a  dimensionless  load  parameter  m  as  follows: 


Rl 

pV 


(10) 


where  R  =  pL  is  the  electrical  resistance  of  the  generator  per  unit 
area.  The  electromotive  force  induced  by  the  generator  is  given  as 
follows: 


<P  =  S(T(0)  -  Tc), 

while  j  is  given  as  follows: 

0  S(T( 0)  -  Tc) 

J  R  +  Rl  R(1  +  m)  ’ 
from  which  we  get 

zr^ 

f  1  +mv  w  ’ 

Therefore,  the  power  output  per  unit  area  is  given  by: 


(H) 

(12) 


(13) 


W=j2R i  = 


mS2(T( 0)  -  Tc)2 

R(1  +  m)2 


(14) 


The  heat  generated  by  the  imbedded  heat  source  in  the  thermoelec¬ 
tric  element  per  unit  area  can  be  evaluated  as  follows: 


(15) 


from  which  rj  can  be  approximately  obtained  as: 
mp2  mZT  mZT 

n  =  ^rt»— w  =  — - —2-  (19) 

Zl  4(  4(1  +m  +  ZTc) 

Eq.  (19)  attains  its  maximum  at  m  =  1  +  ZTC.  The  approximate  max¬ 
imum  conversion  efficiency  is  given  as  follows: 


n°  ~  16(1  +  ZTC)'  (2°) 

As  expected,  Eq.  (20)  shows  that  the  conversion  efficiency  increases 
as  the  thermoelectric  performance  and  heat  source  strength 
increase. 

On  the  other  hand,  when  ZT*  »  1 , 


u  |  zr  (1  +  2m)  1_yzr(l+2m) 

and 

;  yzr(i+2m)  /  zr 

^  ~  1  +  2m  (  y  1  + 2m 

Therefore, 

_  mp2  m  ZT*  m  1 
^“ZT^ZT  l+2m~  1  +  2m  ^  2 5 
which  gives  the  upper  limit  of  the  conversion  efficiency  as 
r\ o  <  1/2  =  50%.  It  should  be  noted  that  this  number  is  not  some¬ 
thing  that  can  be  boasted  of,  although  it  looks  high  at  first  glance. 
In  order  to  obtain  this  value,  the  highest  temperature  in  the  system 
must  be  unlimited,  and  the  thermoelectric  performance  must  be 
infinitely  high.  For  traditional  thermoelectric  generators,  Eq.  (1) 
would  lead  to  the  maximum  efficiency  of  100%  in  such  a  case 
(Th/Tc  and  ZT  — >  (x>).  Eq.  (23)  rather  shows  that  the  present  design 
inherently  has  a  handicap  compared  to  the  regular  design.  The  rea¬ 
son  for  such  loss  of  efficiency  can  be  attributed  to  the  uniform  dis¬ 
tribution  of  the  heat  source.  Thermoelectric  elements,  whose 
temperature  ranges  between  the  highest  temperature  and  the  low¬ 
est  temperature  in  the  system,  generates  heat  in  the  present  design. 
Thus,  the  temperature  of  its  heat  source  cannot  be  characterized  by 
the  highest  temperature  in  the  system,  and  hence  the  thermody¬ 
namic  efficiency  suffers.1  On  top  of  that,  this  number  (50%)  overes¬ 
timates  the  potential,  and  realistic  performance  is  much  worse,  as 
presented  in  Section  3. 


(21) 

(22) 

(23) 


and  the  conversion  efficiency  can  be  simply  given  as: 


3.  Examples 


W 

Qh 


2  (0(0)  -  ecy  -  = mv. 


Solving  Eqs.  (8)  and  (13)  for  p,  we  get: 


I* 


C 

1  +2m 


(16) 


(17) 


where  f  =  1  +  m  +  ZT0C  =  1  +  m  +  ZTC.  In  principle,  substitution  of 
Eq.  (17)  into  Eq.  (16)  leads  to  an  expression  of  q  as  a  function  of  m, 
which  can  be  maximized  to  yield  the  maximum  conversion  effi¬ 
ciency  for  given  geometry  and  material  (k,  q,  S,  p,  and  L  specified). 
The  resulting  expression  is,  however,  too  complicated  to  provide 
any  meaningful  interpretation.  Instead,  we  study  two  limiting  cases 
where  ZT  is  either  very  small  (ZT*  «  1)  or  very  large  (ZT*  »  1). 

When  ZT  -c  1 ,  Eq.  ( 1 7 )  can  be  expanded  as  a  series  in  ZT*  up  to 
the  lowest  order  as  follows: 


In  this  section,  the  performance  of  the  proposed  thermoelectric 
generator  is  analyzed  with  varying  design  parameters.  We  consider 
SrTi03  as  the  material  for  our  thermoelectric  element,  whose  prop¬ 
erties  are  summarized  in  Table  1.  The  values  of  S,  p,  and  k  roughly 
correspond  to  the  values  reported  in  [7],  The  maximum  value  of  q 
is  based  on  the  average  power  density  of  the  90SrTiO3  pellets  with 
7%  excess  Ti02  used  for  the  SNAP-7A  generator  [12],  Note  that  the 
value  of  q  can  be  tuned  down  by  lowering  the  content  of  90Sr  in  the 
material. 

The  performance  of  the  generator  is  evaluated  as  we  vary 
L,  m,  q,  and  Tc.  With  the  necessary  parameters  specified,  one  can 
first  calculate  p  using  Eq.  (17),  obtain  0(0)  using  Eq.  (8),  and  finally 
evaluate  W  and  q  using  Eqs.  (14)  and  (16). 

While  typical  elements  in  conventional  thermoelectric  modules 
are  only  several  millimeters  long,  a  length  scale  of  several 


(18) 


1  Theoretically,  it  is  possible  to  overcome  the  handicap  of  low  efficiency  by 
concentrating  the  heat  source  near  the  hot  end,  which  may  be  a  topic  for  a  future 
investigation. 
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L  =  1  cm;  q  =  0.92  W/cm3 


crosses  represent  the  actual  performance  of  the  generator  and  the  approximate 
value,  represented  by  Eq.  (20),  for  each  of  three  values  of  Tc,  which  are  300  K 
(highest  performance),  400  K,  and  500  K  (lowest  performance),  respectively. 


centimeters  is  more  reasonable  in  our  case,  where  the  thermoelec¬ 
tric  element  also  acts  as  a  heat  source.  In  order  to  generate  a  signif¬ 
icant  amount  of  heat,  a  heat  source  must  possess  a  large  enough 
volume.  For  example,  the  diameters  of  the  90SrTiO3  pellets  of  the 
SNAP-7A  generator  were  reported  to  be  3.4-3. 8  cm  [12],  Cases 
where  1^10  cm  are  considered  in  this  study.  Tc  is  varied  between 
300  K  (room  temperature)  and  500  K  (elevated  temperature).2 

The  performance  graphs  of  typical  cases,  where  q  =  0.92  W/cm3, 
are  provided  in  Fig.  2.  For  a  given  geometry  (I),  material  (S,  p,  k, 
and  q),  and  boundary  condition  (Tc),  the  performance  graph  exhibits 
one  optimal  operating  condition  specified  by  a  particular  value  of  m. 


2  Although  the  present  design  eliminates  the  thermal  interface  at  the  hot  side,  it 
still  needs  to  maintain  a  good  thermal  interface  at  the  cold  side.  The  condition 
Tc  =  500  K  is  used  to  address  the  effect  of  the  thermal  resistance  at  the  cold  side.  As 
shown  later  in  this  section,  a  typical  power  output  is  about  2000  W/m2,  and, 
considering  the  efficiency  in  the  order  of  1%,  a  typical  thermal  load  that  must  be 
removed  at  the  heat  sink  is  at  most  about  0.2  MW/m2.  This  value  is  much  smaller 
than  the  critical  heat  flux  in  the  boiling  heat  transfer  of  water  at  1  atm  (more  than 
1  MW/m2),  which  occurs  at  the  excess  temperature  of  30  °C  or  equivalently  at  the 
surface  temperature  of  130  °C  [15],  500  K  is  significantly  higher  than  130  °C,  and 
hence  the  effect  of  a  finite  heat  transfer  coefficient  through  the  heat  sink  is  partially 
addressed. 


The  operating  condition  for  maximum  power  output  and  that  for 
maximum  efficiency  are  identical,  since  Qh  is  independent  of  m,  un¬ 
like  conventional  thermoelectric  generators  where  the  amount  of 
the  heat  transfer  from  the  hot-side  thermal  reservoir  depends  on 
the  operating  condition,  i.e„  the  load  resistance.  The  maximum  effi¬ 
ciency  of  the  generator  is  very  low  (much  less  than  1%)  especially 
when  ZT  is  small,  as  shown  in  Fig,  2(a).  Increase  in  ZT*  by  increasing 
L  improves  the  efficiency,  as  expected  from  Eq.  (20),  but  Fig.  2(b) 
shows  that  the  actual  value  still  remains  relatively  low  (»1  %)  even 
with  L  =  5  cm. 

The  power  density  is,  on  the  other  hand,  much  better.  With 
L  =  5  cm,  the  maximum  power  density  may  reach  a  level  beyond 
400  W/m2,  which  is  the  power  density  of  photovoltaic  power  gen¬ 
eration  that  can  be  achieved  by  today’s  best  experimental  designs 
(multi-junction  concentrators  with  efficiency  of  about  40%)  during 
the  peak  insolation  hours  [16],  It  should  also  be  noted  that  the 
approximate  value  of  the  maximum  efficiency  computed  by  Eq. 
(20)  is  close  to  the  actual  value,  which  suggests  that  Eq.  (20)  can 
be  used  to  roughly  estimate  the  performance  of  the  proposed  de¬ 
sign,  though  the  difference  between  the  approximate  value  and 
the  actual  value  increases  as  ZT*  increases. 

More  insight  on  the  way  how  the  generator  works  can  be 
brought  in  by  looking  at  a  typical  temperature  profile  inside  the 
thermoelectric  element,  as  shown  in  Fig.  3.  There  are  a  few  things 
to  remark.  First,  the  temperature  profile  is  parabolic,  mainly  due  to 
the  effect  of  the  heat  source  embedded  in  the  thermoelectric 
element.  It  is  a  very  different  profile  from  the  nearly  linear  shape 
typically  observed  in  traditional  thermoelectric  generators  [13], 
Second,  the  maximum  temperature  (Tmra)  occurs  not  at  x  =  0  but 
at  a  point  strictly  inside  (x  =  0.3  cm  in  this  case),  though  the  differ¬ 
ence  between  Tmax  and  T(0)  is  hardly  recognizable  (only  about  1  K). 
This  is  due  to  the  competition  between  the  heat  conduction  and 
the  Peltier  heat  at  x  =  0.  As  one  can  see  from  Eq.  (3),  the  incoming 
Peltier  heat  must  be  compensated  by  the  outgoing  heat  conduction 
in  order  to  maintain  the  adiabatic  boundary  condition.  Since  the 
electromotive  force  given  in  Eq.  (1 1 )  is  dependent  only  on  the  tem¬ 
perature  values  at  the  edges,  the  parabolic  shape  and  the  occur¬ 
rence  of  Tmax  at  a  point  inside  may  have  a  somewhat  detrimental 
effect  to  the  performance  of  the  generator.  A  way  of  alleviating 
such  a  negative  impact  can  be  a  topic  that  should  be  pursued  in 
a  future  study. 

In  order  to  investigate  the  maximum  efficiency  and  power 
density  achievable,  the  performance  indices  have  been  computed 
while  the  design  parameters  (L,  q,  Tc)  are  changed.  For  each  set  of 
design  parameters,  m  is  optimized  to  achieve  the  maximum 


Fig.  3.  A  typical  temperature  profile  inside  the  thermoelectric  element,  q  = 
0.92  W/cm3,  I  =  5  cm,  Tc  =  500  K,  and  m  =  1  +ZTC. 
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Wo  [W/m2];  Tc  =  300  K 
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Fig.  4.  Contour  plots  of  performance  variation  in  cases  where  q  <;  0.92 
W/cm3,  L  .<  1 0  cm,  and  TC  =  300K,  while  m  is  optimized  for  the  maximum 
performance,  (a)  Maximum  efficiency  and  (b)  maximum  power  density.  Curves 
represent  the  contour  lines  of  T(0)  (dashed:  800  K;  dotted:  1000  K;  solid:  1200  K). 


performance.  The  result  for  Tc  =  300  K  and  that  for  Tc  =  500  K  are 
presented  in  Figs.  4  and  5,  respectively.  Obviously,  the  performance 
indices  improve  as  q  and  L  increase,  for  ZT  is  enhanced  in  the 
course.  However,  unchecked  increases  of  q  and  of  L  are  impossible 
for  several  practical  reasons.  First,  as  mentioned  earlier  in  this  sec¬ 
tion,  q  can  only  be  tuned  down  from  its  maximum  value,  since  the 
amount  of  Sr  atoms  that  can  be  replaced  with  90Sr  atoms  is  limited. 
Another  rather  more  significant  constraint  comes  from  the  consid¬ 
eration  of  T(0).  The  increases  in  q  and  L  naturally  result  in  an  in¬ 
crease  in  the  temperature  of  the  inner  core,  which  can  be 
detrimental  to  structural  and  material  integrity  of  the  system. 
Especially,  excessive  temperature  of  the  inner  electrode  may  criti¬ 
cally  affect  the  stability  of  the  system.  Although  noble  metals,  like 
platinum,  rhodium,  palladium  or  iridium,  can,  in  principle,  be  used 
to  produce  high-temperature-stable  electrodes,  there  are  both 
technical  and  economic  challenges  in  the  application  of  these  noble 
metals  [17],  Thus,  the  maximum  attainable  value  of  T(0)  should  be 
considered  as  a  limiting  constraint  here.  For  this  reason,  the 


Fig.  5.  Contour  plots  of  performance  variation.  The  same  convention  as  in  Fig.  4  is 
employed  except  for  Tc,  which  is  set  as  Tc  =  500  K. 


contour  lines  of  three  values  of  T(0)  (800  K,  1000  K,  and  1200  K) 
are  also  given  in  Figs.  4  and  5.  Table  2  summarizes  the  results  from 
Figs.  4  and  5.  The  results  show  that  the  maximum  efficiency  is 
more  or  less  limited  to  the  level  of  2%  or  3%,  while  the  maximum 
power  density  may  reach  around  3000  W/m2. 

As  the  last  example,  we  consider  a  situation  close  to  the 
condition  of  an  actual  radioisotope  thermoelectric  generator.  The 
comparison  is  made  against  the  general-purpose  heat  source  radio¬ 
isotope  thermoelectric  generator  (GPHS-RTG),  which  was  used  as 
the  Galileo  and  Ulysses  power  sources  [9,10],  The  GPHS-RTGs  were 
built  to  deliver  approximately  300  W  of  electrical  power  with  a 
nominal  fuel  loading  yielding  about  4400  W  of  thermal  energy 
[9,10],  yielding  an  efficiency  of  6.82%.  Since  the  thermoelectric 
unicouples  in  the  GPHS-RTG  operated  at  a  nominal  1273  K/573  K 
hot  junction/cold  junction  in  a  vacuum  environment  [9],  we  set 
Tc  =  573  K  and  T(0)  s;  1273  K.  The  expected  performance  of  our 
sandwich-type  thermoelectric  generator  in  such  a  condition  is 
reported  in  Table  3. 

The  maximum  efficiency  of  the  proposed  design  in  the  GPHS- 
RTG-like  condition  is  2.76%,  which  is  lower  than  that  of  the 
GPHS-RTGs.  This  result  is  not  surprising.  As  already  mentioned 
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Table  3 

Maximum  efficiency  and  power  density  in  the  GPHS-RTG-like  condition. 
q  ^  0.92  W/cm3. 
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573 
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8 

2.76 
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earlier  in  Section  2,  it  is  expected  that  our  sandwich-type  thermo¬ 
electric  generator  has  a  lower  thermodynamic  efficiency  than  that 
of  traditional  ones.  However,  since  the  concept  of  embedding  heat 
source  in  thermoelectric  elements  is  new  and  underdeveloped, 
there  is  a  room  for  further  design  optimization,  which  may  yield 
future  improvement  in  performance.  Additionally,  we  note  that 
the  highest  temperature  in  the  entire  generator  in  our  case  is  more 
or  less  limited  to  1273  K,  while  the  highest  temperature  occurring 
in  the  GPHS-RTG,  which  would  be  the  temperature  at  the  center  of 
the  heat  source,  must  be  much  higher.  Accordingly,  thermal  and 
structural  issues  may  become  less  severe  in  our  design,  which 
can  be  seen  as  an  advantage  over  the  traditional  system. 


4.  Summary  and  discussion 

In  this  paper,  a  new  design  of  thermoelectric  generators  in 
which  a  heat  source  is  combined  into  thermoelectric  elements 
has  been  proposed.  By  placing  the  heat  source  inside  the  thermo¬ 
electric  elements,  the  limitation  due  to  heat  transfer  is  only  located 
on  the  side  toward  the  heat  sink,  simplifying  the  overall  heat 
management.  The  use  of  radioactive  isotopes  in  conventional 
thermoelectric  materials  is  proposed  as  a  method  of  creating  such 
elements.  The  inner  electrode  is  sandwiched  between  identical 
thermoelectric  elements,  which  naturally  causes  the  inner  core  to 
act  as  the  hot  side. 

Our  analysis  in  Section  2  has  shown  that  the  performance  indi¬ 
ces,  that  is,  conversion  efficiency  and  power  density,  increase  as 
the  heat  density  of  the  material  is  increased  and  the  thermoelectric 
performance  is  improved.  The  theoretical  maximum  efficiency  is 
shown  to  be  50%.  Unfortunately,  realistic  performance  is  much 
worse.  As  shown  in  Section  3,  the  efficiency  in  realistic  cases  would 
be  about  3%  at  best.  On  the  other  hand,  the  power  density  of  the 
proposed  design  exhibits  a  much  more  reasonable  value  around 
3000  W/m2. 

Although  its  efficiency  is  low,  the  simplicity  of  the  proposed  de¬ 
sign  combined  with  its  reasonable  power  density  may  find  some, 
albeit  limited,  potential  applications,  like  power  generation  for 
control  units  in  isolated  storage  sites  of  radioactive  nuclear  wastes. 
In  order  to  realize  such  potential,  however,  more  serious  study 
beyond  this  preliminary  one  must  be  performed.  Several  technical 
issues  involving  the  proposed  design,  which  have  not  been  ad¬ 
dressed  in  this  preliminary  study,  may  emerge.  For  instance,  po¬ 
tential  radiation  damages  to  thermoelectric  and  other  functional 
materials  may  place  severe  limitation  on  the  heat  density  of  the 
heat  source  that  can  be  embedded  in  thermoelectric  materials. 
Radiation-induced  segregation  and/or  void  formation  [18]  may 
eventually  lead  to  poor  thermoelectric  performance.  Resolution 
of  such  issues  needs  further  future  investigation. 
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Appendix  A.  A  Remark  in  Eq.  (23) 


At  a  glance,  Eq.  (23)  seems  to  lead  to  a  contradictory  conclusion 
that  the  maximum  efficiency  is  achieved  when  m —>  oo,  where  the 
power  output  must  become  zero.  This  is,  however,  only  an  appar¬ 
ent  contradiction.  What  we  have  done  during  the  analysis  leading 
to  Eq.  (23)  is  asymptotically  taking  the  hypothetical  limit  of  infi¬ 
nitely  good  thermoelectric  performance,  in  order  to  identify  the 
absolute  maximum  efficiency  that  can  be  achieved  by  the  present 
design.  For  any  finite  ZT*,  the  maximum  efficiency  is  achieved  not 
at  m  =  oo  but  at  a  finite  value  of  m.  This  can  be  shown  by  expand¬ 
ing  Eq.  (21)  up  to  the  next  order,  as  follows: 


ZT(1  +  2m)  \/ZT(\  +  2m)  , 

C2  “  " 


e 

v'zr(l+2m)J’ 


which  gives  the  following  expression  for  jl\ 

/  ZT*  L  c  j 

^  i/ZT(l  +  2m)J ' 

Therefore, 


~l+2m[  yzr/l  +2m)J 


(26) 


where  f  =  1  +  m  +  ZTC.  Under  the  condition  of  Tc/T  — >  0  and 
ZT*  — >  oo,  the  value  of  m  that  maximizes  Eq.  (26)  is  given  as  follows: 


(27) 


Therefore,  with  a  large  but  finite  ZT*  given,  the  value  of  m  that 
yields  the  maximum  efficiency  also  remains,  albeit  large,  finite. 

The  seemingly  contradictory  phenomenon  is  physically  caused 
by  the  competition  between  the  heat  conduction  and  the  Peltier 
heat  at  x  =  0,  as  shown  in  Eq.  (3).  As  we  increase  j,  the  associated 
Peltier  heat  becomes  large,  making  T( 0)  lower  than  Tmax.  Since 
the  power  output  and  the  efficiency  are  all  affected  by  T( 0),  as  in 
Eqs.  (14)  and  (16),  the  increase  of  the  Peltier  heat  atx  =  Oadversely 
affects  the  performance  indices.  At  a  very  high  ZT*,  where  such  an 
adverse  effect  can  become  significant,  m  must  become  relatively 
large  to  avoid  a  negative  consequence. 
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